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Abstract 

We calculate meson correlators in the e-regime within partially quenched chiral 
perturbation theory. The valence quark masses and sea quark masses can be chosen 
arbitrary and all non-degenerate. Taking some of the sea quark masses to infinity, 
one obtains a smooth connection among the theories with different number of flavors, 
as well as the quenched theory. These results can be directly compared with lattice 
QCD simulations. 



1 Introduction 



In the low energy limit, the dynamics of QCD is described by the pion fields 
which appear as pseudo-Nambu-Goldstone bosons accompanying the spon- 
taneous breaking of chiral symmetry. Chiral dynamics and chiral perturba- 
tion theory (ChPT), play an essential role in understanding the interactions 
among the pions themselves, as well as their couplings with the other (heavier) 
hadrons and sources. 



The fundamental parameters of ChPT are unknown coupling constants in the 
effective theory. But they can be determined by non-perturbative and first- 
principle calculations in the underlying theory, QCD. The most important low- 
energy constants are the chiral condensate E, and the pion decay constant, F, 
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both at leading order. Numerical simulations of lattice QCD give currently the 
most promising approach to achieving this non-perturbative determination. 
Because such computer simulations necessarily are restricted to finite volumes, 
it is important to understand finite- volume effects in the effective field theory. 

Near the chiral limit, the finite- volume effects become increasingly significant 
due to the diverging correlation length of the Goldstone bosons. In particular, 
when the pion correlation length, or the inverse of the pion mass m^r overcomes 
the size of the box L, 

-^«L«— , (1.1) 
Aqcd rriT, 

where Aqcd is the QCD scale, the zero-momentum mode has to be treated 
non-perturbatively and the ChPT has to be performed in a way that achieves 
this in a systematic fashion: the so-called e-expansion [1,2,3,4]. New counting 
rules are needed to order the perturbative expansion in this case. In units of 
the ultraviolet cut-off, the expansion parameter can be defined as e, where, 
unusually, the pion mass m.,^ is not treated as being of the same order as pion 
momentum p. Instead, 

~ ~ e', (1.2) 

while the inverse of space-time volume V and quark mass rriq are being of 
C(e^). The particularly important combination rnqTV is thus treated as of 
order unity. With this expansion, the precise analytical predictions for physical 
observables in the low-energy sector of QCD at finite volume V are expressed in 
terms of the low-energy constants at infinite volume. By comparing numerical 
results at finite volume with these predictions, one can thus extract the infinite- 
volume constants directly from finite-volume simulations, without the need for 
extrapolations of data to infinite volume. The closer one gets to the chiral limit, 
the bigger is the advantage of this approach. 

A few years ago, the predictions for correlation functions of ChPT were ex- 
tended to the cases of both quenched QCD and full QCD at sectors of fixed 
gauge-field topology [5,6]. For the chiral condensate, the studies were also ex- 
tended to the partially quenched cases [7,8]. In quenched QCD these analytical 
predictions suffer in the e-regime from quenched finite- volume logarithms [9]. 
Strictly speaking, such logarithms prevent taking the infinite- volume limit, 
and basically invalidate the whole chiral expansion in this regime for the 
quenched theory. The hope is that in finite ranges of volume, the resulting 
predictions may still have a certain range of validity. For both quenched and 
unquenched cases these computations were restricted to the case of degenerate, 
light, quarks. 

As mentioned above, the great advantage of the predictions for correlators in 
the e-regime is that they appear almost tailored for numerical lattice compu- 
tations near the chiral limit. Indeed, exploratory studies of these correlation 
functions have already shown the great potential [10,11,12,13,14,15,16]. In par- 
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ticular, the possibility of using to one's advantage the role played by fixing 
topology in finite volume has been clearly demonstrated. Also the analytic 
handle one has on the quark mass dependence due to the finite size effects 
in the e-regime has proved helpful in reducing the systematic errors of lattice 
simulations. 

In this paper, we generalize these analytical computations of correlation func- 
tions in the e-regime to the partially quenched theory with separate valence 
and sea quarks, both of which are taken to be non-degenerate. The chiral 
condensate and pseudo-scalar and scalar meson correlators are calculated as 
functions of non-degenerate quark masses, topological charge, and the volume 
of the Euclidean space-time V. In a separate forthcoming publication [17] we 
will present the analogous results for the vector and axial vector channels. 

In practice, the e-regime is not trivially reached in numerical simulations. It 
is therefore important to be able to go close to the chiral limit, but still 
only marginally in the e-regime, while valence quark masses are taken to 
that regime. In fact, although our aim in this paper is the e-regime, our par- 
tially quenched chiral perturbation theory (PQChPT) [18,19] in the e-regime 
smoothly connects all the theories with a different number of flavors as a 
function of the sea quark masses. In this way, our calculation interpolates be- 
tween the e-regime and the more conventional p-regime, and in one kinematical 
regime also mixes the two expansions. For this reason we provide expressions 
not just with the (simpler) e-expansion propagators, but the more general ex- 
pressions. If one of the sea quark masses is taken to inflnity in the A^/-flavor 
theory, it converges to the {Nf — l)-flavor theory. Even the quenched theory 
can be obtained by carefully introducing the flavor singlet fleld before taking 
all the sea quark masses to infinity. Of coiuse, the low-energy constants in 
addition have an inherent flavor dependence that is beyond our control. 

Our results have wide applicability to unquenched lattice QCD studies near 
the chiral limit [20,21,22]. One can choose various valence quark masses with 
a flxed sea quark mass. The partially quenched condensate and meson corre- 
lators can be compared with simulations of heavier quarks which are perhaps 
only marginally in the e-regime (or beyond) while the valence quark mass is 
still in the e-regime. 

This paper is organized as follows. In Sec. 2, we describe the leading contribu- 
tion of the partition function of PQChPT in the e-expansion. We discuss, in 
particular, the exact non-perturbative integral of the zero-modes [23,24] which 
plays a crucial role in deriving both the chiral condensate and meson correla- 
tors in this extended theory. As one fundamental building block of this work, 
the chiral condensate and its 1-loop level correction are obtained in Sec. 3. 
The exact zero-mode integrals in the replica limit are calculated in Sec. 4. We 
also derive a non-trivial identity which follows from the unitarity of the group 
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integrals. In Sec. 5, our main results on meson correlators are presented. We 
plot Nf = 2 connected pseudo-scalar and scalar correlators as examples. The 
conclusions are given in Sec. 6. 



2 The partition function of PQChPT 



Our starting point is the {Nf + A^)-fiavor chiral Lagrangian, 




(2.1) 



where S and F are the chiral condensate and the pion decay constant at 
infinite volume, both in the chiral limit. In the mass matrix 



we have in mind a situation in which the valence quark mass is always taken in 
the e-regime of m^'EV ~ while the physical sea quark mass rrii may vary 

more freely. Unlike standard chiral perturbation theory, U{x) is an element of 
the U{Nf + A^) group, and the flavor-singlet field, ^o{x) = ^^^Tr In [/(a;) , is 
introduced explicitly as a physical degree of freedom with additional constants, 
mo and a. The number of colors is denoted by Nc- As is well known, in this 
partially quenched theory one can normally take the mo oo limit without 
difficulty. In terms of first replicated and then quenched valence quarks one 
is then going from U{Nf + N) to SU{Nf) in a smooth way. Then $o can be 
decoupled from the theory. Of course, trouble arises again if we consider the 
theory in a regime where the sea quark masses have effectively decoupled. 
We will discuss this issue below. 

Separating the zero-mode, Uq, and the non-zero modes, i{x), 



we consider three types of expansion of the partition function in a sector of 
fixed topological charge v. 

(1) Both of the valence and sea quarks are in the e-regime: 



M = diag(mi,m2, ■ ■ ■ ,mNf,m^, ■ ■ ■ ,my), 



(2.2) 



U(x) = C/oexp(iv^^(x)/F) , 



(2.3) 



4 



'Nf+N 



[niy, {rrii}) 



U{Nf+N) 



dUod^ det Uq exp 



Ti[M^Uo + MUl] 



+ J d'x i^-^TT[d,^d,^]^ + 



(2.4) 



where the mo — >■ oo hmit is taken and the singlet $o is decoupled from 
the theory. 

(2) The sea quarks are in the p-regime (but light enough that we can still 
disregard effects of the singlet field) : 



U{Nf+N) 



dUod^ det Uq exp 



+ / d'x {-\TT[d,id,i] - (-|) Tr[A^^^] + 



(2.5) 



Here the valence sector is expanded as in e-expansion, while the sea sector 
is considered in the p-iegime. Note that the mass term in the valence 
sector, nivT,^'^ / F"^ , is of O(e^) but not ignored here, in order to see a 
smooth transition to the p-regime. 
(3) The sea quarks are heavy: 



U{Nf+N) 



dU^d^ det Uq exp 



+ j d'x ^-}-Tr[d,^d,i] - 



2 

-) 

FV 



MUl] 



2N, 



(TrO^ 



a 



(5.TrO^ + 



(2.6) 



As the sea quark mass increases, new terms (the non-zero mode's mass term, 
and the singlet fields) come in. Thus, Eq.(2.6) is the most general form. 

In the following, we implicitly restrict ourselves to sectors of small enough 
fixed topology v for the expansion to be valid [25]. Indeed, this is an essential 
assumption so that in Eq.(2.5) and Eq.(2.6) one needs only the mass term 
Tr[7W^2] instead of Ti\M{Uq + Uq)(^]/2\}\. 

For the zero-mode integrals, one needs exact formulae of the group integrals 
over U{Nf + A^) and means of taking the replica limit. As described in detail 
in ref. [5], if we wish to consider correlation functions with external valence 
quarks we must embed the valence quarks in a theory with A^ replicated 
quarks in total (of which A^ — A^ do not couple to the external sources), and 
then take the limit A^ 0. Alternatively, one can consider a theory with A^^, 
additional bosonic flavors of common mass m^. It is easy to understand in the 
quark determinant picture that this limit is equivalent to the replica limit: 



We thank P. Hernandez for stressing this point. 
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Hin det(i:' + m)^^ det{D + + Ji) j det(i:' + m^) 



^ ^ det(D + m^)^" ^ ^ 

Here the left hand side is the rephca prescription, while the right hand side is 
the prescription of the graded formalism. 

The zero-mode partition function of n bosons and m fermions are analytically 
known [23,24], 



det[nl ^X+i-i(/^i)]j,i=i,- 



■n+m 



(2. 



where /Xj = miTiV . Here ^'s are defined as Jj,^j-i{iXi) = {—iy~^Ki,+j_i{ni) 
for i = 1,- ■ - n and J'u+j-i{fJ,i) = Ii,+j-i{lJ>i) ior i — n + 1, ■ ■ ■ n + m, where 
Ii, and are the modified Bessel functions. Of particular importance is the 
case (n, m) = (1, Nf + 1): 



( K,{x) h{y) h{z^) h{z2) 

i(^i) ^2-^1^+1(^2) 

X'^K^+2{.X) y^I^+2{y) zfl^+2{Zi) Zp^+2{Z2) 



X det 



(2.9) 
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where x — mfjUV {nib denotes the bosonic quark mass), y = m^'EV and 
Zi — rriiEV. One notes that 

^^^2^1,l+Nfi^\y, = 2^0,N,iU^}), (2.10) 

and therefore, 

- I^^y9x^i,i+Nfix\y, {zi}) = lim ^^.^^^^^^^^(a^ly, {zi})- (2.11) 
It is also remarkable that 



^l,l+Nf{x\yj {Zl, Z2, - ■ ■ , Zj-i, Zj — » 00, Zj+i, . . . , ZNf}) — 

^l,l+{Nf-l){x\yj {Zl, Z2, - ■ ■ , Zj^i, Zjj^.1, . . . , Z^f}) 



X 



Zj) 



(2.12) 
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which is consistent with the intuitive notion that the Nf — 1-flavor theory 

can be obtained in the hmit of large sea-quark mass, rrij oo, up to a 
normahzation factor Zq^^'^-^^^zj) . This decouphng is of course 

more generaL For example, in the case of {zi — > oo (for all i)} one obtains 
Zi^i, the leading partition function of quenched chiral perturbation theory. 

For different valence quarks (with non-degenerate valence masses) , we will also 
need Z2,2+Nf- 



22,2+Nf{XuX2\yi,y2,Ui}) = 

1 

4 - xDivl - yf) nS(^.? - yl){zf - yl) TluiMl ' 4 
^ K,(xi) K,{x2) h{yi) IM iu{zi) ^ 



X det 



-xiK^+i{xi) -X2K^+i{x2) yil„+i{yi) y2h+i{y2) ^i-^i/+i(^i) ••• 
xlK^.^2{xi) xlK„+2{x2) yjlu+2{yi) yllu+2{y2) 4+2(^1) ••• 

V ••• ••• ) 



(2.13) 

Let us now define the propagator [26] of the fluctuation field ^: 

Pmki){x -y) = {^^jixKkliy)) 

5ikSkiA{Mi,\x -y)- SmG{MI Ml^\x - y) {i ^ j) ' 

(2.14) 

where the indices can be taken both in the valence and sea sectors. 

Here M^. = (m^ + mj)J:/F^ and 



1 pipx 

^(MiM-yZ^^, (2.16) 

G(M^ b) = i y ^"^K + V)/^c .2 16) 

^(/).1 + E ^i.\y- - (2.17) 
/=i P +^^-^ff 
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Note that for small sea-quark masses, mo can be taken infinity, 



1 pipx 

G{Ml M^,\x) = - E V ■ (2-18) 



Conversely, in the quenched limit of taking the sea quark masses to infinity, 
G becomes 

G(M„M,,W = -X.p-;^^p-^. (2.19) 



3 The chiral condensate 



At tree level, the partially quenched chiral condensate is obtained by the 
logarithmic x-derivative of the zero-mode partition function Eq.(2.9) followed 
by the y ^ x limit. 



d 

-1 



I 



X det 



I„{x) I^{zi) h{z2) 



-d^{xK^+l{x)) Xl„+i{x) Zilj,+i{zi) Z2lu+l{z2) 

d^{x'^K^+2{x)) X'^ly+2{,X) zjl„+2{Zi) Z^h+2{Z2) 



\ 



(3.1) 



where x = niyT^VjZi = miJ^V, and we use Eq.(2.10) for the denominator 
and the minus sign in the first line is due to the derivative with respect to 
the bosonic flavor. The quantity T,^^{x, {zi})/E has two properties that follow 
immediately from the corresponding statements about the partition functions: 

(1) When one of the sea quark masses is large, the partially quenched chiral 
condensate reduces to that of the {Nf — 1) -flavor theory; 
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Nf flavors 

})/S. (3.2) 

^ V ' 

(Nf-l) flavors 

As a particular case, one recovers the fully quenched condensate [5] when 
all sea quark masses are large, i.e. when — > oo (See Fig. 1). 
(2) If the valence quark mass m^, is equal to one of sea quark masses nij, it 
is equivalent to the j-th flavor quark condensate in the full theory: 

lim^. 5?!(|M) ^ lnZ]^^({..}) ^ ^'^"'"''^({^^}) (for any j). 

(3.3) 

The first property follows directly from Eq.(2.12) and the fact that 

^or^^V^o,! - 1 (3.4) 

in the large mass limit. The second property can be shown explicitly by using 
Eq.(2.10) and noting that the product HiJii^i ~ x"^) Hkyji^l ~ ^"j) always 
antisymmetric under a swap oi x ^ Zj. Together with the equation 



/ 



det 



K^{x) d^.Iu{,Zj) 
-xK^+i{x) dz.{zjh+i{zj)) 
x'^K^+2{x) d,j{z]lu+2{Zj)) 



j-th column 

LiZj) 
ZjIv+\{Zj) 
ZjI„+2(Zj) 



\ 



I 



det 



Kyix) I^{zj) 
-xKi,+i{x) Zjly+i{zj) 

X^K^+2{x) Zp^+2{Zj) 



/ 

J-th column 

dz^h{zj) 

d,^{zjl^+i{zj)) 
d,.{z]I,+2{zj)) 



\ 



(3.5) 



which holds for any j, the statement follows. 

For the normalization factor of the correlation functions to be computed in the 
following sections we need to evaluate the partition function to the given order 
in the expansion. This essentially boils down to an evaluation of the one-loop 
correction to the partially quenched chiral condensate. This correction can be 
calculated in a standard manner: 
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SPQ,l-loop(^^ {^,}) 1 Id 

E = TV i^^o Nd^^ ^iv,+^K, {m.}) 



/ 1 JV/+iv ^ 



1 ^ 

Xl^(E[^o + ^o']..,)^;^°°^ (3.6) 

i V ^ 

where (■ ■ • )^ denotes the integral over the ^ field, while the zero-mode integral 
(■■■)c/o^°°^ can be calculated as above but due to the vacuum bubble, the 
arguments x and Zi are shifted, at 1-loop, as 



m^Yy = x^x''^ and m{i:V = Zi ^ zf , (3.7) 
where 



= iim 1 

X N-.0 F^NV 



^= hm I 1 - 



^ Nf+N N \ 

2^ T.{J d'x^if^f^)A . (3.8) 



x'^ 




1 




= 1 - 




X 




F2 


^eff 




1 




= 1 - 






F2 



With the Fcynmann rules given in the previous section, the ^ integral can be 
performed, and we obtain [27] 



'Nf 



Y.K{Ml^\Q)-G{MlMm \ ■ (3.9) 



The corrections are UV divergent, and thus need regularization. Note that 
each of x and Zi receive different one-loop correction in general. In practice, 
however, the following three special cases are of our interest: 

(1) Both of the valence quarks and the sea quarks are in the e-regime: 

One can to this order set — for all i and j, and take the mo oo 
limit, which leads to 

z^^ 1 A^? - 1 - 

— = ^ = 1 - j-,^^Hm (for all i) . (3.10) 
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This correction is equivalent to that of A^/-flavor full theory. In dimen- 
sional regularization, 



(3.11) 



is obtained where Pi is known as the shape coefficient [4]. It depends only 
on the shape of the 4-dimensional Euclidean space-time volume. 
(2) The sea quarks are in the p-regime while the valence quarks are still in 
the e-regime (the m^'s are heavy but much smaller than the QCD scale, 
^QCd): 

One can take the mo oo limit but should keep m^'s finite, which leads 
to 





= 1 - 


1 






X 


F2 


^eff 




1 




= 1 - 








F2 



Nj 



Nf 



\ 



Ea(m2|o)--E 



1 



(3.12) 



Note that a double pole contribution appears in x'^^ /x, as an effect of the 
partially quenching. 
(3) All the sea quark masses are heavy, rrii ^ Aqcd: 

In this case, mo cannot be large, but one can take the m^ — > oo limit, 
which leads to 

which agrees with the quenched result. 
To summarize this section, the chiral condensate to one-loop order is given by 

™cfT 

5.PQ,i-ioop^^^ {;,,}) = EPQ(a;^^ {zf}) — , (3.14) 

where the analytical functional form of T,^^{x, {zi}) is given by Eq.(3.1). When 
all the quarks are in the e-regime, the one-loop correction is, to this order, 
constant, and can simply be taken into account in the Lagrangian by shifting 
S according to the above prescription. The chiral condensate in the infinite 
volume limit, S, and all the other low-energy constants, are of course expected 
to depend on the number of fiavors. Matching conditions [28,29,30,31,32] can 
ensure smooth connections between theories with different number of flavors. 
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N^2, v=0, x=m^Ey=2.0 fixed 

1 I 




1000 



Fig. 1. The sea quark mass dependence of the 2-flavor partially quenched con- 
densate at 1/ = (solid). The valence quark mass is fixed to myTiV = 2.0. As 
expected, we obtain a smooth curve crossing the full 2-flavor theory result (dashed) 
at mi = m2 = m^, and converging to the quenched limit (dotted). 



4 Zero-mode integrals in the partially quenched theory 



4-1 U{Nf + A^) group integrals in the replica limit 



In the e-regime, the integral over the zero-mode, Uq (which for simplicity of 
notation will be denoted by U in this subsection) has to be done exactly. This 
is the central difference between the e-regime and the p-regime, and we are 
fortunately able to perform the required group integrations exactly. Because 
these technical aspects are so important for the calculation that follows, we 
first give a detailed outline of how the integrations have been done, and how 
we employ the replica formalism in this context. As in Ref. [5] the needed 
group integrals are conveniently obtained through the identities that follow 
from the fact that [det{J]VM)]'~''Z'^^,_^_j^{x, {zi}) is a function of M only. 
Since the group integrals are known for diagonal sources, this basically solves 
the problem. We begin, however, with two simpler cases. They do not require 
any special techniques, and can be evaluated straightforwardly by making use 
of the graded partition function: 
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\{{Uvv + Ul)) = \imdy\nZl,^^^{x\y,{z,}) = ^^^^^^^^^^^ (4.1) 



_ d,Y.l%x,{z,}) ASPQ(x,{^,}) 
E E 

where the second term of Eq.(4.2) is defined by 



(4.2) 



(4.3) 



If we put X = Zj (for any j), this simply amounts to removing the fermion 
determinant of quark species j. From the definition (4.2) we then immediately 
get 



d,El<i{z,,{zi}) AEl^{z„{z,}) 



2 



d.,- — ,^^1 , (4.4) 



corresponding to the result in the full theory without partial quenching. Con- 
versely, in the limit 2;^ — > oo (for all i), 



d,El^{x, {z^}) d,T.r{x) 



E {zj}— »oo E 

,,4. (4.5) 

E {zi}-»oo X-^ 

we recover the results of the quenched theory [5] . Here, in a hopefully obvious 
notation, we have denoted the chiral condensates in the full and quenched 
theories by E^,""'^ ^^^'^'^ and E?"^ (See Eq.(3.2) and Eq.(3.3)), respectively. 



We next consider a purely imaginary source iJ (with J real) on a diagonal 
{v,v) element of M., 

T.VM + J = diag(2;i, Z2, ■ ■ -, (a; + iJ), x, ■ ■ ■ , x). (4.6) 
' » ' ^ , ' 

Nf N 

We note that 

Um^ det(EyAl + J) = det(EyA^) x (1 + iJ/x), (4.7) 
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and 



Nf 



N 



where A = v^j?+j2 gy ^j^g chain rule, J-derivative and A-derivative are 
related through 



d_ 
dJ 



Jd_ 



0, 



^2 



aj2 



J=0 



ld_ 
X dX 



(4.9) 



Prom the above equations, we obtain (denote Z^^_^_j^{x, {zi}) by Z for sim- 
plicity). 



^^det(El^A1 + J)-^Z|j=o = det(El^A1)-'^ f -- + ^^^j = 0, 
det(El^A^ + J)-''Z| j=o = det(EyA^)-" ( ^^^^^^ + ^ 
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J=0> 



ll^det(EVA< + J)-'^ZU=, 
Z X aA 



This gives us the useful identities 



V x^ xIj 



(4.10) 



ZdJ^ 



V 

X 



J=0 



xE 



2 ■ 



(4.11) 
(4.12) 



In order to calculate the meson correlators in PQChPT, we also need ma- 
trix elements which have different valence flavor indices. Por example, for the 
disconnected correlators, we need 



= DPQ(a;i,a;2,{^.}). 



(4.13) 



To derive the analytical expressions for these and others closely related, we 
flrst consider two purely imaginary sources Jv^vx — iJ\ and Jv<2,v2. — iJi along 
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the diagonal: 

T,VM + J = diag(2;i, Z2, ■ ■ •, (xi + iJi), (xs + 1.12), x, ■ ■ ■ , x), (4.14) 
^ , ' ^ V ' 

Nf N 

The rephca hmit of this determinant is 

hm^det(EyAt + J) = det(ET/M) x (1 + iJi/xi)(l + iJa/xa), (4.15) 

and {EVM + J^iEVM + J) takes the form 



{EVM + jy{EVM + J)^dieig{zlzl--;XlXl,x^,--- (4-16) 



Nf 



N 



where Aj = yxf + (i = 1,2). Again, J-derivative and A-derivative are re- 
lated through the chain rule. 



Ji=0 



A, 



0. 



Ji=0 



From these equations we get 



(4.17) 



1 

ZdJM 



detiT.VM + jyZ 



= det(EyA^)-'^ 



1 d'Z 



X1X2 Z dJ\dJ2 



0, 



(4.18) 



Next let us consider a real off-diagonal source, J^viv2 = for which the deter- 
minant is unchanged in the replica limit: 

hm det(EVA1 + J) = det{J:VM) . (4.19) 

Now {J:VM + jy{T,VM + J) can be diagonahzed as 

{EVM + jy{EVM + J)^ diag(^^ ■ ■ •, A^, Ai, • • • , x^) , (4.20) 



Nf 



N 



where 



A 



± 



( J2 + xl + xl) ± JJ^ + 2J2(a;f + xl) + (x^ - x^f 



(4.21) 
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If we assume that xi ^ X2 (the special case xi X2 [5] can be taken as a 
hmiting case afterwards, see below), we obtain the following relation between 
J-derivatives and A-derivatives: 



^2 



dp 



1 



d 



+ 



d 



X2 



A+=a;i 



(4.22) 



A_ =X2 / 



We now use 



XlTT. X2-T— I Z 



dX 



A-^ —xi ,A_ —X2 

, SPQ(xi,{^.}) EjQ(x2r[^^ 



1 



(4.23) 



Note that a purely imaginary source ^viv2 — iJ gives the same results. Thus, 
one obtains 

i<(^.» ± UUf) - - . (,.4, 



As a check, if we take the mass-degenerate limit X2 Xi we recover 
hm +r7t ,2,_^^( K'^{^l.{^^}) , d,,Y.l%x,,{z,}) \ 



(4.25) 



which is obtained from the formula in the degenerate case in Eq.(4.22), 



^2 



aj2 



d 



d 



d 



d 



4x1 V^'^+ A \d\+ d\. 



(4.26) 



A±=a;i 



where A± = {^J P + Axl ± J)/2. 

We finally put two real sources on the off-diagonal elements, JviV2 — Ji ^^'^ 
•Jv2Vi = ^2- In the replica limit the determinant becomes 

lim^det(El^M + J) = det(EVA<)(l - JiJs/xiXs), (4.27) 

and {TyM + J)'^{TyM + J) diagonahzes as 

{TyM + J)\TyM + J)^ diag(^2, zl ■ ■ •, A^ , Ai, x^, • • • , x^), (4.28) 



Nf 



N 



where 
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A 



± — 



2 

(4.29) 



The relation between J-derivative and A-derivative can be worked out as 
above. Assuming again that xi ^ X2 (the degenerate case can also here be 
recovered by taking the limit Xi X2 afterwards), we get 
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dJidJ2 



d 



J =0 x{ — x^ \ dX^ 



d 



Xl 



A+=a;i 



9A_ 



A_ =X2 1 



(4.30) 



In the same way as above, we thus find that the two equations 
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det(SyA^ + jy^Z 



dQi{T.VM)- 



V 



+ — 



1 S'Z 



x\X2 Z dJ\dJ2 
[^2^ - a^i^l det(SVA< + J)-'^ZU,=., 



X 



+ 



dX. 
1 



det(SVA^)-'^ 



X\X2 X^ — X2 



. SPQ(xi,{^,}) El^{x2,{zi}y 

2 ^21(^2 Xl 



lead to 



(4.31) 



1 f ^l'^{x„U^}) ^l'^{x2,U^}y 

2 ^2 ^; 2:1 

X-^ X2 \ ^ 



where we used that purely imaginary sources J'viv2 = ''"^1 Jv2Vi = 
the same result. We summarize all pertinent formulas in appendix A. 



(4.32) 
iJ2 give 



J!i..2 The Unitarity Formula 



From the above formulae and the requirement of unitarity one finds in the 
replica limit. 
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N+Nf 



Nf 



N-1 



Nf 



(4.33) 



where we have used that the N rephcated quarks are all degenerate: Xi — 
rriy^'EV — X. One then obtains a non-trivial identity, 



AEPQ(x,{^,}) 



(4.34) 



It is not difficult to show that the left hand side of Eq.(4.34) actually vanishes 
in the quenched limit, cf. Eq.(4.5). 

It is also interesting to consider the limit of the full A^/-flavor theory of 
Eq.(4.34), in which all the valence and sea quarks are degenerate: x — zi — 
Z2 — • • • — z. In that degenerate limit we have 



E^^J_ Er(^^'^)(.i,.2,---) 



Nf^ E 



Z-i=Z2 = ---=Z 



Zl=Z2 = ---=Z 



+ {Nf-l) 



(4.35) 



Zl=Z2 = ---=Z 



(4.36) 



where we have used the fact that the 92^.S[)^"*^*)'s are independent of i or j in 
the generate case. We thus have 
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21=Z2 

/ 



Nf - 1 



Zl=Z2 = 



(4.37) 



Similarly, one obtains 



zi=Z2=---=z 



Nf-ld^El'^ix, {z,z,...}) 



Nf 



+ 



1 l^uHz) 



Nf 

x=z J 



~ Nf[ E 
{N]-l) d,J:l<i{x,{z,z,---}) 



Nf 



dz,Er'^''^''\{z,,z,z,...}) 



(4.38) 



The unitarity equation for the degenerate case {x — zi — Z2 — ■ • ■ — z), 
Eq.(4.34) then finally becomes 



Nfj: 



4.39) 



Zl=Z 



a relation which is useful when simplifying the expressions for the meson corre- 
lators in Scc.5. We finally note that Eq.(4.39) is consistent with an analogous 
formula obtained from Schwinger- Dyson equations [5] in the full theory. 



5 Meson correlators 



In this section we present the detailed analytical predictions of partially quenched 
scalar and pseudoscalar correlation functions in the e-regime of QCD. Our 
calculation is done to the lowest non-trivial order in the e-expansion, and by 
taking limits we can check that we recover both the fully quenched and Nf = 2 
results at sectors of fixed topological index i/ as reported in ref. [5]. 
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5.1 Partially Quenched Correlators with Two Valence Quarks 

Of most immediate interest are the correlation functions with two hght valence 
quarks, representing the u and d quarks in QCD. The physical case of low- 
energy QCD with two very light quarks (the u and the d quarks) and one 
heavier quark (the s quark) is an example where one can be in a mixture of 
the e-regime (with respect to the two light quarks) and the p-regime (with 
respect to the s quark). Even if the two hght quarks correspond, for the given 
lattice volume, to the e-regime, it may be convenient to recycle the lattice 
configurations by considering light valence quarks that are also in the e-regime, 
but just at different mass values. This requires a comparison with the results 
of PQChPT in the e-regime that we shall present here. 

Furthermore, with the resulting formulae as building blocks one can calculate 
various other different types of meson correlators. They can correspond not 
just to different valence and sea quark masses but also to an arbitrary num- 
ber of valence and sea quark ffavors. For example, the 5'[/(3)-singlet scalar 
correlator is obtained by 

{S'{x)S'm = EE(g.,g.,(a;)g.,g.,(0))d + i:(g.,g.,(x)g,,g,,(0))e, (5.1) 

i j i 

where the expressions in the r.h.s. are defined below. Note that the contri- 
butions inside the sums can differ due to the non-degeneracy of the valence 
quark masses. 

The connected scalar and pseudo-scalar correlators are defined by 

(0))c 

(5.2) 

(fe75fe(a^)fe75fe(0))c 

X {U{0)y,y, - U''{0)y,y, + U{0)y,y, " U''{0)y,y,)) , 

(5.3) 

where Vi denotes the valence flavor index with the quark mass rriy. = Xi/T^V. 
Note that Eq.(5.2) and Eq.(5.3) have the same structure as the iso-triplet 
correlators in the 2-flavor theory of which masses are m„j and ruy^. 
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The disconnected correlators are similarly defined 



(0)). 

(5.4) 

{(iviimv^ {x)qv2i5qv2{o))d 
'(y{x) 

(5.5) 



Let us begin with the disconnected scalar correlators. To O(e^), we find 



)P{12){2l){x) 



2F2 
"2F2 



V2 / eff \ ^ 



E2 



- C/i.J(C/.... - UlJ)G{Q. 0|x), (5.6) 



2F2 

We have consistently set m^^ = = in the NLO contributions, and 



eff eff 1 / ^/ 

^ = ^ = 1 - ^ E A(M^|0) - G(0, 0|0) I , (5.7) 



which of course needs regularization. Note (■ ■ ■ )^ indicates the shift Xi — > 
xl^ and — > zf^ in the arguments of the Bessel functions. 

In the same way we obtain the disconnected pseudo-scalar correlation function 
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J]2 / ^cff \ 2 



2F2 

2F2 



X] 



(5.8) 



For the connected correlators, we get 

(0))c 

•£2 / ^cflF \ 2 



/ cflF \ 2 ^ 

T V^J (^^^^^^ + + ^^^^^^ + ^4.J(^^2.i + ulj) " 



-loop 



■^2 / 



2F2 



E [(^.n^J.i)^(2i)(i2)(^) + {u.2iU!j 



-Jp^iUviViUy2V2 + Ulv,Ul^V2) [^(12)(21)(a:;) +^(21)(12)(2^) 

X;2 

"4P^((^^l^2 ~ t^i^l)(t^^2«l - Ulv2))P{22){ll){x) 

-^{iUv2v, - UlJiU,,,, - UlJ)P^nK22){x), (5.9) 



{Qvi75qv2{x)qv275qvi{0)) 



^ \Xi 



{UviV2 ^V2Vl) {UviV2 ^V2Vl){^V2Vl ^VlV2) 



1— loop 



^2 / 



2F2 
"4F2 

E^ 
"4i?2 

E^ 
"4i?2 



E [(t^.ut^i)^(20{i2)(^) + -P(ii)(a)(2^) 

i 

{Ulv2 + (t^U)')^(22)(22)(x) + {Ul,, + (t/4,j')P(ll)(ll)(x) 

{Uv,viU^2V2 + ^i^i^i,;^) [Ai2)(2i)(a;) +P(2i)(i2)(a;) 



{{Uv2Vl + Ulv2)iUv,V2 + U12V^))P{11){22){X), 



(5.10) 



where the summation over the flavors, denoted by X^^, has to be taken carefully. 
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For example, in the sum 



i 



P{2i){i2) {x) 

i 

+ {Uv,V2Ul^v^)P{22){22){.x) - {U^,v^Ul^v,)P{22'){2'2){x) 
Nf 

+ T.iUv^iU}jPi2i)ii2)(x), (5.11) 
i 

the index with a prime i' is treated as a different ffavor from i-th qiiark but 
with the same value of the quark mass. With this, we find the correlators 

{qv^qv2{x)qv2qvi{0))c 
y^2 / cfr\2 

+ 2^2 (^■"l^'l^ill)! ~ ^V1V[U^,^^^ + UV2V2UI2V2 ~ ^'"■2'"'2^v'2V2 

" i^lJlVi ^V2Vl ^V2VlUl^y2 — {Uv\V2 ~ ^V2Vl)^V2Vl ~ ^iill2 ) ) ^ ^\x) 

"2F2 



+^ E(^^i^< + U,2iUl2)HMl\^)^ (5-12) 



(fe75g^;2(a^)fe75gw(0))c 

4 \ 2^1 / \ ^^''^ ^V2Vl) ^ \'^ViV2 V2Vl)\^V2V\ ^V'iV2) / 

+ ^(t^.l.2t^.U + t^.2.lt^i.. + {Uv,V2 + C/.U)(C/.2.1 + t^i„J)G(0, 0|X) 

Y.{UvnUl, + t/.,.t/;,)A(M2 (5.13) 
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where we have expUcitly kept the sea quark mass, in order to possibly apply 
these equations outside of the the e-regime. 



With aid of the formulae collected in Appendix A and use of the defining 
equation (4.13) for D^^, one obtains 



T? 2 (x^K^^ix^^iz,}) x{Ll%X2,{zi}Y 



F'^xi-xl\ E 



A(0|x) 



{qvilbqv^ (a;)fe75fe (0))<i = — . 

\Xi J Xi X2 

2 (x_^^{x^^ _ x{Ll%x,,{zC^) \ 
F^xl-xl\ E E J^^UIxj 

0:2, {^.})G(0, 0|a;), (5.15) 



■sPQ(xf,{zf}) EPQ«,{;.f})^ 



E^ MEPQ(xi,{z,}) AEPQ(X2,{^,}) 
2F2 I E E XI X2 



+2DPQ(a;i,a;2,{^i}) + 



a;ia;2y 



A(0|a;) 



T? \ 4 ^EP Q(.Ti,{^.}) ^ EPQ(X2,{^.}) ' 



2F2 [xi + X2 V 



G(0,0|a;) 



.2 iV/ 



2F2^ 



2 /xiEPQ(xi,{z,}) z,e1^'^(^^'^)({z,})' 



2 /x2EPQ(a:2,{^,}) z,Er'(^^'^)({z,}) 



- z] 



A{MUx), (5.16) 
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(g^;l75fe(2^)fe75g^>l(0))c = 



cfr\ 2 



X 



1 ( s^ Q«,{^f}) ^ sPQ(xf,{zf}) ' 



cflF I ^eflf 



\ 



2F2 



-2L>PQ(xi,X2,{z,})- 



+ 



2u^ 

X1X2 



+ 



™2 /Y»2 
Xi X2 



A{0\x) 



2F2 



G(0,0|x) 



^2 A'/ 



2F2 



E 



2 /a;iSPQ(a;„{^,}) ^,Er^"^^'^'\{^.}) 



2 /X2EPQ(X2,{^,}) ^,Er^^^'^'\{;.,}) 



- z] 



A{MUx). 



(5.17) 



We have here three types of (^-correlators, 



1 pipx 

^(0W = T7E^> (5-18) 

G(0,0|x) = -|: . (5.20) 

Note that if the sea quarks are much smaller than the cut-off of ChPT, but 
still in the p-regime, one can take the mo 00 limit. It leads to 



1 pipx 

= — A(0|x) + E 1 + (5.21) 

where the well-known double pole contribution appears due to a mismatch of 
the sea and valence quark masses. Further simplification is possible when all 
the sea quarks are in the e-regime, 

A{Miy\x) A{0\x), G{0,0\x) -^A{0\x). (5.22) 

Nf 
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An interesting special case is the degenerate limit xi — X2, where the above 
formulae become 



X 



efr\2 



{qviqvi ix)qv,qvi (0))d = -s 



AE^Q(xf , {zf}) 



Xi 



x) 



S2 2^/2 ^ 



F2 



G(0,0|x), 



(5.23) 



{qvil5qv^ (a;)g^i75g«i (0))^ = S' 



,eff\ 2 



\ Xi J {Xi^ 

F2 I E 



XiE / 



2E2 AEPQ(xi,{^,}) 
F2 E 



G(0,0|x), 



(5.24) 



efr\ 2 



Xi 



a.EPQ«,{^f}) 



^^^^^'(-'^^^>)g(o,o|.) 



F2 



xiE 



-2^(01^) 



F2 X? 



.2 N; 



2F2 



E 



4 (xj:l^{x^,{z.}) z,Y.f'^''^^'\{z,}) 



(5.25) 



(gi,i759i,i (a;)g^i 75^1,1(0)), 



-eM^ 



EPQ(a;f,{zf}) , 2E2AEPQ(a;i, {;.,}) 



V xi / xf E 

2E2 a.EPQ(xi, {;.,}) - 



+ 



F2 



F2 

V2 ^ 



G'(0,0|x) 



2F2 



E 



:riSj:^(.ri.{..}) 



A(0|x) 



A(iW|Ja;), 
(5.26) 
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5.2 Singlet and flavored meson correlators 



To compare with lattice QCD data, the case with arbitrary A^„-flavor valence 
quarks is interesting. In particular, can be different from Nf, the number 
of sea quarks. In this paper we have aimed at the e-regime predictions, and 
of course we always take the valence quarks to be in that regime. With the 
qualification mentioned in section 2, one can also be more general, and consider 
the sea quarks to be just marginally in the e-regime, or even entirely in the 
p-regime. The finite volume is a x T box, where L and T are the spacial 
and temporal extents, respectively. The zero-momentum projections of singlet 
and {N^ — l)-plet correlators are then obtained 



I d^x{S%x)S''{0)) 



E2 



efr\ 2 



d^El^ixf, {zf}) AEPQ«, {zf}) 



2F2 

E^ 
"2F2 

2F2 



X 



b{t/T) 



N E^ 

j 



XiE 



4 [x^J:l^ix„{z,}) z,Ei^'"^''^'^\{z,}) 



a{t/T) 



Xi - z] 



Cj{t/T), 
(5.27) 



j d^x{P\x)P''{0)) 
= - j Ny{q^^'yr,q^^{x)q-u^-f5q^^{0))c + N^{qvi^5qvi{x)qvi^5qvA0))d 



,eff\2 



sPQ«,{^f}) 
,xf S 



2F2 



4iV,ASPQ(a;i,{^,}) 



E 

-2N 



,,d,El^{x„{z,}) El^{x„{z,}y 



2F2 



E xiE 
d^El^ix„{z,}) _ AEPQ(xi, {z,}) 



a{t/T) 



N E2 
2F2 ^ 



4 / iiEP«(xi, {^.}) ^^Er'^'-^'Cl^i}) 



(5.28) 
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Kb 
2 



cff ' 



d^El'^{xf,{zf}) 2E2i/2 



Xi 



F2 xl 



a{t/T) 



2E2 EPQ(xi,{z,}) 



F2 



xiE 



E2 

2F2 ^ 



xiE^Q(xi, {zi}) ZjEu 



full(7V/,i) 



({^.}) 



- Zj 



c,{t/T) 
(5.29) 



y d'x{P\x)P\Q)) = -^-f j d'a;(g„,75g.i(x)g,,75g„,(0)), 



V 



2 

2E^ a.EPQ(a:i,{z,}) 
E2 



EPQ«,{zf}) , 2E2AEPQ(xi,{z,}) 



xfE 



+ 



F2 



h{t/T) 



2F2 



E 



4 /xiEPQ(xi,{z,}) z,i:'f^''^^^\{z,}) 



x\ — z^ 



where a{t/T),b{t/T) and Cj(t/T) are defined by 



a{t/T) 



cAt/T) 
(5.30) 



T 








2" 




2. 


' 12 



a{t/T) = j (i^xA(0|x) 

6(i/r) = y (i3a;(5(o, o|x) = y ^^^^E 



^tpx 



1 T 







]--] 




2j 





'TV 



r ^1/ 24 



Tj \T 



30 



1 



axiAym^^lx) 2M,„ sinh(M,-,r/2) Mf,T 



T 
2 



t 1 



1 

12 



3'V 

T) \T j 30 



+o(mM, 



(5.31) 



(5.32) 



(5.33) 



where M? = m^T./F'^ and Mf. = 2miE/F2. 
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In Fig. 2, Fig. 3 and Fig. 4, we plot, as examples, the flavored pseudo-scalar 
and scalar correlators (we simply denote Eq.(5.29) and Eq.(5.30) as {S"'{t)) 
and (P"(t)) for a = 6 ). We use S =(250 MeV)^, F = 93 MeV, xf/xi = 1 
and L = T = 2 fm as inputs. 



5.3 Ward-Takahashi identities 

In this subsection, we check that the above results satisfy the Ward-Takahashi 
identities under the chiral rotation of the degenerate Ny valence quarks. For 
the singlet chiral rotation, one obtains 

{{d^Alix) - 2mP\x) - 2iN^uj{x)0{Q)) = -{60{0))6{x), (5.34) 

for any operator 0{x), where 



Ny Ny 

v=l v=l 

oo{x) = ^TrF^,F;,(a;), (5.35) 

and 60 denotes the chiral variation of O. Note that Eq.(5.34) holds not only 
in ^-vacuum but also in a fixed topological sector. The identities for 0{x) ~ 
P^{x) and 0{x) — u!{x) and their integration over the volume give an equation 



J d^x(P°(x)P°(0)) 



Ny {S%x)) _ Ny ^N,El'i{xf,{zf})xf 



m?V m rrrV mxi 

(5.36) 

where S^{x) — Yl^li Qv{x)qv{x), which coincides with Eq.(5.28). Note here we 
used 



J dt a{t/T) ^ Jdt h{t/T) ^ j dt Cj{t/T) = 0. (5.37) 

In the same way, the flavored identity 

{{d^A'^^ix) - 2mP''{x))O{0)) = -(TO(0))<5(a;), (5.38) 
for 0{x) = P^ix) and 0{x) = d^A^^ix) gives 



29 



A 
V 



A 



CO 
V 



45 
40 
35 
30 
25 
20 



45 







20 







v=0, m^= 5MeV fixed 



m^=5 MeV 
m^=3 MeV 
m^=l MeV 




0.5 1 

t[fm\ 



v=0, m^= 5MeV fixed 



1.5 



m^=5 MeV 
m^=3 MeV 
my=l MeV 



0.5 



1 

t[fm\ 



1.5 




Fig. 2. Pseudo-scalar (top) and scalar (bottom) flavored correlators with both 
valence and sea quarks in the e-regime, here for 1^ = 0. The sea quark masses are 
fixed to 5 MeV. We use S =(250 MeV)^, F = 93 MeV, xf/xi = 1 and L = T = 2 
fm as inputs. 
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Fig. 3. The same as Fig. 2 but with v = 2. Note that the pseudo-scalar and scalar 
correlation functions almost sum up to zero. 
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Fig. 4. The sea quark mass dependence of pseudo-scalar (top) and scalar (bottom) 
flavored correlators with the valence quark mass fixed at 3 MeV. The parameters 
are chosen the same as for Fig. 2 and Fig. 3. Here we assume that we can still ignore 
the mo and a terms for sea quark masses up to ~ 15 MeV. 
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which is also consistent with Eq.(5.30). Here the operator with the superscript 
"a" has the form of qy{x)jT'^q.f,{x) with some gamma matrix 7, as in the 
conventional notation, where T" denotes the a-th generator of SU {Ny) group. 

One can also confirm the similar identity [33] 

(5.40) 

is consistent with Eq.(5.29). Since the asymptotic form of the partially quenched 
condensate in the chiral limit is known [7] 

Z_i Xi 

it is not difficult to see that the known "quenched" identities for the Ny = 2 
case 



/r n 4 7/ 4-;y 

d'x l{P%x)P\0)) - 4{S^{x)S^m] ~ -^F - (5.42) 

/ d'^x [(P"(x)P'^(0)) - 4.{S\x)S\m (5.43) 
J rriyV 

also hold in the limit my — > in this partially quenched theory. 
5.4 Quenched and full degenerate Nf flavor limits 



In this subsection, we show how to reproduce the known results of both the 
fully quenched theory and the unquenched theory by taking the Zi 00 limit 
and the x = zi = Z2 = ■ ■ ■ = z limits, respectively. The valence quarks are 
chosen to be degenerate in the e-regime. 



First consider the quenched limit of connected correlators Eq. (5.25) and Eq. (5.26) , 



(5'«(x)5»(0))«- = Jim l(g.,?,,(x)^,,?,,(0)), 



cff\ 2 



X 



Xi 



d^E^ixf) 



2E 



2P2 



Xi XiE 



(5.44) 
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(P'^(a;)P'^(0))^- = lim --(g„,75fe {x)q.,l,qv, (0)), 



Noting 



X 



eflfX 2 



2F2 



-2(i + ^U(o|.) + H^5p£l)e(o,ow 



X 



1/ 



(5.45) 



A(0|x) 



1 pW^ _ 1 

G(0,0|a;) = lE 



1 / e^^^^mQ e*^^Q; 



+ 



(5.46) 



in the quenched Umit, one can see that Eq.(5.44) and Eq.(5.45) agree with the 
quenched results in [5]. 

Next we construct = 1 singlet correlation functions in the quenched limit, 



(5°(x)5°(0))''"^= lim {{M,M^,,q,M)c+{MvAx)^vAvM)d) 



cfr\ 2 



X 



Xi 



2F2 



E + + {xff\ 
4i/2 2E«"^(xi) 29^E^""(xi)^ 



x\ 



+ 



a:iE 



t2 



A(0|a;) 
(5.47) 



(P°(a;)P°(0))«"^ 

= lim ((g^,l75g^l(a;)g„l759m(0))c - (^^;i75?w (a^)^^;i75?^;i (0))d) 



V 



efr\2 /5]<?«e(^efr) 



(Xf )2 



2P2 



E 

4i/2 2Ef^(xi) 2a^E«"^(xi)\ , , 
4 + — ^4r^ + " ^ ^ A(0|x) 



+ 



'4a,Er(^i) 



+ 4 + 



,TiE 

4z/2\ - 

— G(0,0|x) 



(5.48) 



which are also equivalent to the results of [5] . 

The full degenerate A^^-fiavor limit [x ^ z = Z\ 
correlators in the e-regime are 



^2 



) of the connected 
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z . 



2E 



2Z.2 2s;y^(z) /^S^QMMi. 
^ + lv^"^H ^ — 



}) 



A(0|a:), 



(5.49) 



(P-^lxjP^lO))^"" = 75^.1 75^.1 (0)) 



c|a;=zi=Z2"-=z 



z / 2^<=f^E 

E2 r 2AEPQ(^,{z,z,-.-}) 29.EPQ(a;,{^,z,---}) 



^ 2F2 



a,EPQ(x,{z,^,---}) 



+ 



zE 



A(0|a;). (5.50) 



Note here that we set Ml = M^^ = 0, G'(0,0|a;) = A(0|a;)/iV/, i:7{z)/ll is 
the full degenerate A^/-flavor condensate defined by Eq.(4.35), and z'^^ is given 
by Eq.(3.10). To ehminate the partially quenched expression, E^*^/E, we use 
Eq.(4.34), Eq.(4.38) and Eq.(4.39) to obtain 



(5"(a;)5"(0))^''" = 



1 

'Iff 



2{N] -l)\z 

E;'^(^^f^) 



Nf 



effj] 



+ 1 + 



2(iV| - 1)F2 



1 



^E 
A(0|x), 



(5.51) 
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V2 



2{Nf - 
+{Nf + 2) 



(iV|-4)^ + 3iV,^ + iVj-4 



V 



which agree with those in the full A^/-flavor theory. 

In the same way, one can see that the singlet correlators 



A(0|x), (5.52) 



,eff 



-R 



,AEPQ(^'^ff,{^^f^,^-ff,---}) 



Y? 
2F2 



4(iV2 - 1) 



E^(^) 
^E 



A(0|x), 



(5.53) 



(P°(x)P°(0))-^"" = -A^;(g,, 75^.1 (a;)g.i75g.i(0))c-A^/(g.i 75^.1 (2:)g.i 75^.1 (0))d 



2F2 



4(iV2 



a.EPQ(x,{z,z,---}) 

^"^^ E 



A(0|a;), 



(5.54) 



are also consistent with the known expressions in [5]. 



6 Conclusions 



In this paper, we have discussed partially quenched chiral perturbation theory 
(PQChPT) in the e-regime, and in the mixed e and p-regime. 

Using the 1-loop improved chiral condensate and its derivative as building 
blocks, we have calculated various zero-mode group integrals in the replica 
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limit. These integrals are necessary for the computation of mcsonic correlation 
functions in the partially quenched theory. We have also derived a non-trivial 
identity which is a consequence of unitarity of the graded or replicated group. 

With these zero-mode integrals and the Feynman rules for the non-zero modes, 
we have calculated the mesonic correlation functions for both connected and 
disconnected pseudo-scalar and scalar channels with non-degenerate quark 
masses, both of the valence and sea kind. Among others, our results can be 
applied to the mesons that consist of two non-degenerate valence quarks. For 
a demonstration, we have plotted the flavored pseudo-scalar and scalar cor- 
relators with a realistic choice of input parameters. As expected, they show 
a non-trivial valence (sea) quark mass dependence with a fixed sea (valence) 
quark mass. 

These meson correlators were shown to have the correct quenched and degen- 
erate full iVj-flavor theory limits. We have not addressed the implicit flavor 
dependence of S, or F, the fundamental parameters in the infinite volume 
limit. In order to complete the smooth connection among the theories with 
different number of flavors, one has to match the value of them [28,29,30,31,32] 
The flavor dependence is expected to be rather weak, but the matching is in- 
teresting and important for the future works. 

Our results are useful for the analysis of unquenched lattice QCD simulations 
in many ways. The various valence quark masses can be used for each set 
of fixed sea quark masses. Even if the physical pions are just barely in the 
e-regime, one can put the valence pions very safely in the e-regime and com- 
pare numerical data with our formulae for the partially quenched correlation 
functions. 
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A Summary of group integrals 



Here we summarize the group integrals in the rephca hmit, which are necessary 
for the meson correlators. See Sec.4.1 for the details. The formulae for one 
valence index are, 



\m. + ui))=^^^^^, (A.i) 

\m.-Ul))^~ (A.3) 



hm _u^f) = -^M^^-l .A 4) 



2 



For two valence indices, 



+ UlJiU,,,, + UlJ)=Dl^{x,,X2, {z,}), (A.6) 

Im,., - t/i„J(t/.... - UlJ) = (A.7) 

4 X1X2 

{U.,v,U,,.,) + {Ul,,UlJ=2Dl'^{x,,X2, {Zi}) + , (A.8) 

X1X2 



where D^^ is defined in Eq.(4.13). Similarly, 



2 {^V\V2^V2V^ 2 ^^'"^VlU^^y^) 

±1 /^^^ SPQ(a;i,{z.}) ^^ EPQ(x2,{z.}) ^ 



Xi X2 



(A.9) 

^R%. + (C/U)') = 0, (AlO) 
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^{{UviV2 i ^i;2fi ) (^■U2fl ^ ^iiD2)) 

1 / El^{x,,{z,}) 




] 



(A.11) 



E 



are obtained. 
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